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The tight-binding model of quantum particles on a honeycomb lattice is investigated in the pres- 
ence of homogeneous magnetic field. Provided the magnetic flux per unit hexagon is rational of the 
elementary flux, the one-particle Hamiltonian is expressed in terms of the generators of the quantum 
group Uq{sl2). Employing the functional representation of the quantum group Uqi^sh) the Harper 
equation is rewritten as a systems of two coupled functional equations in the complex plane. For the 
special values of quasi-momentum the entangled system admits solutions in terms of polynomials. 
The system is shown to exhibit certain symmetry allowing to resolve the entanglement, and basic 
single equation determining the eigenvalues and eigenstates (polynomials) is obtained. Equations 
specifying locations of the roots of polynomials in the complex plane are found. Employing numer- 
ical analysis the roots of polynomials corresponding to different eigenstates are solved out and the 
diagrams exhibiting the ordered structure of one-particle eigenstates are depicted. 



I. INTRODUCTION 



The problem of quantum particles in the presence of a 
periodic potential and a uniform magnetic field has been 
the subject of intensive studies for decades. Azbel [l| was 
the first who pointed out that the spectral properties of 
two-dimensional lattice particles have a sensitive depen- 
dence on the flux through the plaquette. This observa- 
tion has been exploited later by Hofstadtcr [2| who found 
the exotic structure of the one-particle energy spectrum 
of planar particles on square lattice in magnetic field. 
The same study was extended later for triangular lattice 
Q , generalized square lattices 0, Q and for the honey- 
comb lattice [y, 0] ■ These studies firmly established the 
fractal structure of the aforementioned energy spectrum, 
whose rich and complex nature originates from the pres- 
ence of two, not necessarily commensurate periods. The 
first is given by the lattice structure and the second is de- 
termined by the magnetic field. The relevant parameter 
which determines the band structure is the ratio $/<I>o 
where $ is the magnetic flux per elementary plaquette, 
and $0 — 27r(?i/e) is the magnetic flux quantum. 

In 1994 Wiegmann and Zabrodin pointed out [g] that 
the Hamiltonian responsible for the original result of Hof- 
stadter is closely related to the quantum group Uq{sl2) 
(for mathematical treatment see [^]). Namely, it was 
shown, that the Hamiltonian is expressible in terms of 
X^ generators of quantum group Uq{sl2)'. '3-C = X^+X^ 
with a deformation parameter q determined by the ap- 
plied magnetic field. Employing the functional repre- 
sentation of Uq{sl2) in the space of polynomials, the 
Harper equation was reformulated into the functional 
form where the one-particle eigenstates appeared as poly- 
nomials. The zeros of polynomials unambiguously spec- 
ifying the one-particle wave functions were shown to be 
determined by the Bethe ansatz equations Q- The fact 
that eigenstates are related to certain polynomials associ- 
ated with the quantum group Uq{sl2) and with the Bethe 



ansatz equations provides with the possibility to system- 
atically study the structure of eigenstates. It should be 
noted that so far the studies were mainly concentrated 
on the structure of spectrum. 



Hatsugai et al. 
tions derived in 



10| investigated the Bethe ansatz equa- 
and found out that the zeros of those 
polynomial are not scattered randomly, but are located 
along concentric circles on a complex plane, thereby ex- 
hibiting nontrivial and ordered structure of eigenstates. 

Alongside with the eigenvalue problem on a square lat- 
tice, the analogous problem has been considered for a 
honeycomb lattice as well [y, \J\ . The common observa- 
tion is that the energy spectrum diff^ers from the one of 
square lattice and is equally highly nontrivial. After the 
connection between the quantum group Uq{sl2) and the 
problem of square lattice particles in magnetic field was 
found out in |8|, the similar approach has been applied 
to honeycomb lattice by Kohmoto and Sedrakyan |ll| . 

In contrast to the square lattice, the honeycomb is not 
a Bravais lattice, but consists of two interpenetrating tri- 
angular lattices ("A" and "_B" sublattices) with one lat- 
tice point of each type per unit cell. Therefore the one- 
particle Hamiltonian possesses additional (2 x 2)-matrix 
structure. In terms of the generators of Uq{sl2) it appears 
as 



'K = 







X- 



X- 







(1) 



and carries the anti-diagonal structure with respect to 
(2 X 2) -matrix indices. This causes certain inconvenience 
(eigenvalue problem consists of two entangled equations) 
which can be overcome by considering 'K'^ instead of "K 
since the former is diagonal in (2 x 2) -matrix indices. Tak- 
ing advantage of this observation, the authors of [Il| ob- 
tained the analogue of the Bethe ansatz equation derived 
for square lattice in Q. The Hamiltonian (1) describes 
particle hoppings from a site to the three nearest neigh- 
boring ones. The squared operator "K^ contains second- 



order terms describing double-hoppings A ^f B ^f A 
and B —)' A —)' B. As a result, the Bethe-like equations 
derived in [ll| look quite cumbersome, and the nesting 
procedure has been applied with the aim to reduce them 
to the set of simpler equations. 

In this paper we develop the distinct approach to the 
eigenvalue equation for the Hamiltonian (1) by revealing 
its novel symmetry. This symmetry allows to carry out 
certain analytic calculations right for 'K with no necessity 
of involving J{^. Based on this symmetry we reduce the 
system of entangled equations to a single equation thus 
avoiding the necessity of involving !H^. The obtained 
single equation corresponds to the elementary hoppings 
A <^ B and appears perceptibly simpler than the one 
derived in [ll| from J{^. 

The paper is organized as follows. In the forthcom- 
ing section we comment on the derivation of one-particle 
Hamiltonian from the tight-binding model on honey- 
comb lattice in magnetic field. In Section III we rewrite 
this Hamiltonian in terms of the quantum group Uq{sl2) 
and discuss the representations of Uq{sl2)- Section IV 
deals with Harper equation in polynomial representation. 
Since the honeycomb consists of two triangle Bravais sub- 
lattices, the Harper equation appears in the form of the 
system of two entangled equations. This system is shown 
to possess certain type of symmetry which permits to de- 
couple them and to obtain a single equation describing 
the elementary process of nearest neighbouring hopping. 
Solutions to this equation arc polynomials with real co- 
efficients and Bethe-like equations determining the zeros 
of polynomials are derived. In Section V the results of 
numeric analysis are presented and the locations of zeros 
in the complex plane are depicted exhibiting nontrivial 
and ordered structure of eigenstates. 



II. QUANTUM PARTICLE ON A HONEYCOMB 
LATTICE IN MAGNETIC FIELD 



Honeycomb lattice is built up of two-dimensional array 
of hexagonal unit cells of the side a/v3, with atoms at 
the vertices. Such a structure is encountered in solid state 
physics in various crystals, while the ideal realization of 
the two-dimensional honeycomb lattice is graphene |12| . 
The unit cell is a rhombus of the side a with angles 7r/3 
and 27r/3 at its vertices. Each unit cell consists of two 
atoms, A and B as indicated in figure 1. 

We study the tight-binding model on a honeycomb lat- 
tice with the nearest neighbouring hoppings only. In the 
presence of a homogeneous magnetic field the Hamilto- 
nian under consideration is given by 



H 



Z^ 



,-nn.ir) J 



{r + dn)cA{r) 



h.c 



(2) 



where the sum with respect to r is implied over the sites 
r = jiai + J2a2. Here cl^{r) ic^{r)) and c^g{r + 5„) 
{cg{r + Sn)) are the particle creation (annihilation) op- 
erators on site r of the sublatticc " A" and on site r + 6n 
of the sublattice " B" , respectively. 

Magnetic field B is included in the Hamiltonian via the 
Peierls phases 



ln{r) 



Adl, 



(3) 



where the vector-potential is taken in the Landau gauge 
A={-By,0). 



J 





FIG. 1: Honeycomb lattice (left) consists of two Bravais sublattices A and B. The A-sites are located at r = jiai +J2a2 where 
11,2 = ^(±1, V^)o; and ji,2 are integers. The three B-sites, nearest to a given A-site are located at r -I- ^1,2,3. Brillouin zone 
(right) is arranged in the rectangular form. The vectors 01,2 set by aidj = 2n5ij determine the structure of the Brillouin zone. 



Rewriting (2) in the Fourier representation we find 



H = / \cUk)cAk) 



-thai t 



lc}5(fe)c^(fc + fco) 



-ika2 „t 



c_g(fc + A;o)c^(^) dk ^ h.c 



(4) 



where CA{k) and CB(fc) are the Fourier transforms of 
CA{r) and Cg{r + (5„) respectively, and the integration 
covers the first Brillouin zone (FBZ). 

The vector fco is related to the magnetic field 



27r $ , , 



(5) 



where $ and $o a-re the magnetic flux per elementary 
hexagon and the magnetic flux quantum, respectively. 
We consider 



i^ ~ TV' 



(6) 



where v and N arc coprimc integers, and concentrate on 
the odd values of N. 

Then the Hamiltonian (4) can be presented as follows 
(see Appendix A) 



/mbz 



J{(fc) 



X+{k) 



x-{k) ^ 







(7a) 



(7b) 



where ^(fc) is a (2A^)-component column, and the inte- 
gration covers the magnetic Brillouin zone, which is the 
A^'th part of the first Brillouin zone. 

The N X N matrices X^{k) are given by 

X+{k) = e-*'="i^1'Q + e-»'=°2Q;3, 



+ikai rit « I „+ika2 otnt 



X-{k) == e+^^^iQT^ + e+'^^^P 



(8a) 



(8b) 



where 



' f • 


•00- 


1- 


• 


0- 


• 


0- 


• f 


.10 0- 


•0 0. 



(9) 



and 



with 



Q = diag(g\g2,...,g 



+iT!{v/N) 



g = e 



III. THE QUANTUM GROUP f/,(s/2 



(10) 
(11) 



Using the relation Q^ /3Q = q^Q/3Q^ we find 

[X+,X-]=i'{q-q-^){K~K~^), (12a) 



KX'^K-^ = q^^X 



±2 v± 



(12b) 



where 



K{k) = qe 



-i-ik{ai—a2) 



QPQ^P, 



-Vfc^ = «-ip-*'=("i-"2)otr)«tr)t 



K-\k) = q 



/3TQ/3T 



(13a) 



(13b) 



The relations (12) constitute the definition of the quan- 
tum group Uq{sl2) which is the deformation of usual sl2 
(see [l3| ) with q referred to as the deformation parameter 
(we adopt the normalization of AT* which slightly differs 
from the standard one). Note that particular values of k 
and ai_2 do not affect the structure (12). 



A. Cyclic and highest weight representations 

For every complex q ^ ±1 the quantum group Uq{sl2) 
possesses the highest weight representations analogous to 
those of ordinary sl2- Furthermore, when q is the root of 
unity, the so called cyclic representation also does exist. 



Due to q 



N 



(— I)'' the cases when v is either 



even or odd slightly differ from each other. Hereafter we 
discuss the even values of u. The case of odd values leads 
to identical outcome, and the corresponding technical de- 
tails are collected in Appendix B. 

Introduce the states ijji,'4>2, ■ ■ ■ , ipN hi the form of the 
A^-component columns 

i^,=N-i{q\q'^,q'\...,q^^y. (14) 

These states form the orthonormal complete set 



TV 

E 

n=l 



(V'J)«(V',-)" = k 



N 



(^j)m(l/'i)n = Smn, 



(15a) 



(15b) 



and can be used as a basis in the space of representation. 
The action of matrices /3 and Q on the states (14) is 
expressed by the following relations 



/3-0- = q^^. 



PHj = 9 ■''/'j> 



Qi'j = i'. 



j+ii 



)tl/,- =-07-1, 



(16a) 



(16b) 



(16c) 



(16d) 



where the cyclic identification ipN+j = V'j is understood 
in (16c) and (16d). This is possible due to q^ = 1. 



Using (16) wc find 

X+i^j = q~--e-i'^''y%+^^j+i, (17a) 



X-^j = 5+5e+5^'=«"tjV^j„i, 



K^, 



K-'i^j 



_ „+ikxa 



^2J- 



q--'i)j 



e-'^^'^q-'^^ipj, 



where 



tj =e+3'^'-'^QJ-5 



g— 2fcia„— j+ 2 



(17b) 



(17c) 



(17d) 



(18) 



Due to (17a) and (17b) the operators X^ and X~ can 
be regarded as rising and lowering ones respectively. 
From (17a) and (17b) we find 



DetX± = qT^AfgT^wy3fc„a(^^^^ . . . ^^y 



(19) 

For those values of k^ when tj ^ for any j, the 
operators X possess only nonvanishing eigenvalues. In 
other words, there is neither highest nor lowest weight 
states. Hence (8) and (13) form the cyclic representation 
of (12). Acting on -0j by X+ we obtain i^j+i, and so on 
until wc arrive to ipM- The subsequent action takes us 
back to i}}!- The similar is true for X~ . 

On the other hand, for the value of k^ set by e+^'^="° = 
±iq2-io TffQ have tj^ = and the cyclic representation 
turns into the highest weight one, e.g. setting £5*-'=" — 
iq~2 we find ti = tN+i = 0. Then (17a) and (17b) yield 
X~^ipN = and X^ipi = so that f/'jv and -01 appear as 
the highest and the lowest weight vectors respectively. 

The important comment is in order: manipulating with 
quasi-momentum, the cyclic and the highest weight rep- 
resentations can be transformed one into another. Such 
an interplay has been known for a long time [IJj. 



Then the relations (20) can be written as 

X+f{z) = e-i''^''q-izf{q-^z) + 



+ o n^'ajQ, -j- -Ty 



+ e'^2'^-'°'q+2zf(qz), 



(22a) 



X-f{z) 



„+ik^a TT 1 



q 2z 'J{qz)+ 



rfcc.o„+i^-l f(„-l 



+ e-2'^-'^g+2z-i/(5-iz), (22b) 



Kf{z) = e+--^/(q^z) 



K-^!{z) = e-^'^-^fiq-^z) 



(22c) 



(22d) 



These relations determine the functional representation 
of Uq{sl2). Functional representation for the odd values 
of V is constructed in Appendix B. 



IV. THE EIGENVALUE EQUATION 

Wc study the eigenvalue equation 

l + X-{k) ' 

l + X+{k) 





T 


= E 


' ^ ■ 




.C. 




.c. 



(23) 



Employing the functional representation (22) the 
eigenvalue equation (23) can be rewritten in the form 
of two coupled equations 



B. Functional representation 

Consider an A^-component vector / = (/i, /2, • ■ • , Jn)- 
Writing out the action of X^, K, K~^ on / in the com- 
ponent form we obtain 

{X+f)„ = e-''"''q"-\f„-i + e-"=°^g"/„+i, (20a) 



^(z) + e+^''-°g+5 2^(5z)-h 



e-^''-''q-2z^{q-'z) ^ EC{z), 



C(z) + e+^'=-°g-^z-iC(g^) + 



(24a) 



'ifcao „ — n+1 . 



iX-f)„ = e+'^°ig-"/„+i + e+"^"^g-"+\/„-i, (20b) 



~ 9 fcj;fl „~r"9 



z-\{q-'z) = Eaz). 



(24b) 



iKf)n = e 



_ +ife(ai-02 



-1 f\ _ -ifc(ai-a2 



/n+2, 



\fn- 



(20c) 



(20d) 



where the identification fN+n = .fn is assumed. 

Introduce the interpolating function of the complex 
variable f{z) such that 



/„ = /(e-*^'^''"g"-^). 



(21) 



If any pair (^, (^) is the eigenvector corresponding to the 
eigenvalue E = X, then the pair (^, —() is the eigenvector 
corresponding to the eigenvalue E = —X. 

The main disadvantage of (24) is the entanglement of 
^(z) and C{z). This difficulty has been overcome in [11| 
by "squaring up" the equations (23): repeated applica- 
tion of :K to (23) gives out {I + X-){I + X+)£, = E^^ 
and (I + X+){I + X-)C = E^C which in the functional 
representation appear as 



i^-E'n^) 



-i-ik^a 



qiio^z) + e+*''"°(9+22; + q~2z~^)£^{qz)+ 



+ e-*'^-°(q-'z + q+3z-l)e(g"^z) + er'^''''q-^i{q-^ z) = 0, 



(25a) 



(3-i?2)C(. 



ik^a —1 



C(?'^) 



e+i'=-°(g+22 + 5 2z 1 



)C(92)+ 



1 "!j; 



'{q-2z + q+2z-^)C{q-^z) + e-'''''''qC{q-^z) = 



(25b) 



Due to the "square up" trick some portion of infor- 
mation encoded in (24) is lost in (25). Here wc propose 
essentially different approach allowing to avoid this draw- 
back (disadvantage) of the " square up" procedure. The 
key point of our consideration is the observation, that 
the system (24) is invariant under the following transfor- 
mation 



then find 



X+/,(z)=i,+i/,+i(z), 



X-f,{z)^t,f,.,iz), 



Kf.iz] 



- «+2j 



./.(^), 



az) ^ i^zra~z-') 



c(z) ^ i^zra-z-') 



(26a) 



(26b) 



K-'Mz) 



-''fd^), 



(29a) 



(29b) 



(29c) 



(29d) 



where the parameter ui is set by g"e"*'''^^° = — 1. This 
symmetry allows to express the solutions to (24) as 



where tj has been defined by (18). 

Setting 6+5'=-" = iq^ we find X'/o = X+fN-i = 
yielding that the subspace of (A^— l)'th order polynomials 
f{z) = cq + ciz + - ■ • + cjv_iz is the invariant subspace 
of the group Uq{sl2). Then the equations (24) appear as 

C(z) + igz^qz) - iq'^zi{q-^z) = EQ{z), (30a) 



' az) ' 



f{z) 

[ ±{tzrf{-z-^) J 



(27) 



C(z) + iz^Hiqz) - iz-\{q-^z) = ES,{z), 



(30b) 



where f{z) satisfies the equation 



and admit the solutions in the form of [N — l)'tli order 
polynomials. These polynomials are real, i.e. with real 
coefficients. This statement can be verified by expanding 
^(z) and C^{z) in a series of z^ and rewriting (30) for the 
expansion coefficients. 

The invariance (26) holds for (30) provided cu ~ 2J = 
N — 1 while the equation (28) turns into 

+ e-i'^^^'q-hfiq-'z) = A(zz)"/(-z-l). (28) f{z) + iqzj{qz) - iq-'zf{q-'z) ^ {^lY Xz^^ f{~z-'), 



^+ifexa„+i 



f{z) + e+2''^^q+2zf{qz) + 



The signs "±" in (27) correspond to £' = ±A, respec- 
tively. 

We have thus reduced the system of equations (24) 
to a single equation (28) while £,{z) and C{z) are now 
expressed via the unique function f{z). In Appendix C 
we show that the system of equations (25) appears as 
a corollary of (28) meaning that the later is the basic 
equation relevant to the problem under consideration. 

We proceed to study the equation (28) at special val- 
ues of momenta reducing the functional representation to 
the polynomial one. Consider the monomials fj{z) = z^ 
forming the basis in the space of analytic functions. We 



(31) 
generating N eigenvalues Ai, . . . , Xn- The factor of (— 1)"^ 
can be absorbed in A and therefore we may drop it in 
what follows. 

The equation (31) is derived for v = even. For v = odd 
we come to the similar functional equation, but with the 

i7r(iV-i^) 

deformation parameter q — e " instead of q (see 
Appendix B). Therefore, provided N = odd and ly = odd 
we have N — v = even, i.e. q^ = 1. For that reason 
below we discuss only ly = even. 

Since f{z) is a polynomial we can write down 



2J 



f{z) = Y[{z~z,), 



(32) 



j=i 



where Zi, Z2, • ■ • , Z2J are the zeros of f{z). 

Substituting (32) into (31) and taking z 
obtain 



tZr, 



2.7 , 2J _i 

n l + qZnZj ^ -p-j- L + q ZnZj 



I ZnZj 



which is the honeycomb analog of the Bethe ansatz equa- 
tion obtained in [8| for square lattice. Due to the fact that 
(30) is soluble in terms of real polynomials, the solutions 
to the equation (31) are also real polynomials. Corre- 
spondingly, if z„ is a root then z* must be the root as 
well. This comes in agreement with the invariance of (33) 
under the complex conjugation. 
Rewriting (31) as 



2J 



2J 



A^ni 



■ 1 - + ^^j 



iqz 



n 



qz — Zj 



2J _^ 

iz T-r q z — Zj 



-I- ZZ.i 






-I- ZZo 



and taking the limit z — > we find 



(34) 



2J 



A=n 



j=i 



2J , 



j=i 



2," — iq + - 
9. 



2,7 

n ^^- + o(^ 



where the right hand side must be z-independent. Then 

2,7 



A=n 



(35) 
rhen 

(36) 



j=i 



and 



V (zj + — ) =ig+— . 



j=i 



(37) 



The analogous " sum rule" does exist in the case of the 
square lattice. Namely, the functional relation derived in 
H leads to 



2J , 



j=i 



= 0, 



(38) 



where Zj are the zeros of the square lattice polynomials. 
Taking now the limit z — >■ cx3 in (34) we find 



A^ 



2 7 

l + *(<z-,-i)^z, + ofi). 



i=i 



(39) 



where again the right hand side is z-independent yielding 

2J 



\^^\-rl{q-q-^)Y.. 



(40) 



i=i 



Using (37) we can rewrite (40) also as 



3-?^ 



2,7 



(41) 



We have thus derived the relation among an eigenvalue 
and the roots of the corresponding polynomial. The set 
of roots is determined by (33). Once the roots of a poly- 
nomial are known, the corresponding eigenvalue can be 
calculated using any of (36), (40), (41). 

Remark that the aforementioned invariance of (30) 
may be interpreted in terms of special conformal trans- 
formation 



z ^ w = w(z) 



1 



(42) 



Under this conformal map a quasi-primary field ^(z) 
transforms as 



$(z) -> $'(w) 



dw 
dz 



-h 



$(z) 



(43) 



where h is the corresponding conformal weight [15| . 

Taking h = —J and expressing non-primed fields via 
the primed ones we find 



$(z) =z2'^$'(w), 
$(gz) = q-^z^-'-^'iq-^w), 



$(g- 



qz''-'^ (qw). 



(44a) 
(44b) 
(44c) 



Using (44) and rewriting the systems (30) in terms of 
primed fields we obtain 

C{w) + iqwCiqw) - iq^^wC {q^^w) = E^'{w), (45a) 



C{w) + iw-^^'{qw) - iw-'^^'{q~^w) = EC'iw). (45b) 

Hence, we may identify ('{z) = ^(z) and ^'(z) = C(^)i 
i.e. C{z) = z2"'^(-z-i) and ^(z) = z^''C{-z-'^). 



V. NUMERICAL ANALYSIS 

We discuss the results of numerical analysis on how 
the roots are distributed over the complex plane. These 
studies are helpful for visualization of the peculiarities 
hidden in the Bethe-like equation (33). 

We first present the case of ^ = ^ where wc have 89 
eigenvalues and the corresponding 89 polynomials deter- 
mined by (31). Each polynomial is of order of iV—1 = 88, 
i.e. there are 88 roots for each eigenvalue. For demon- 
strative purposes we present several characteristic cases, 
exhibiting how the root locations in the complex plane 
depend on eigenvalues. For the eigenvalue of maximal 
magnitude the roots are arranged along two spirals re- 
lated via the complex conjugation (figure 2a). For the 
eigenvalue of less magnitude some of the roots (out of 
total 88) regroup into the new branch of circular form 
(figures 2b and 2c) . For the eigenvalue of minimal mag- 
nitude the spirals convert and the roots become arranged 
along two ovals (figure 2d). 
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FIG. 2: Root distributions in the complex plane for — = gg- There are 89 eigenvalues and the 89 corresponding polynomials, 
each of order of A'^ — 1 = 88. For each eigenvalue we have separate set comprising 88 roots of the corresponding polynomial. 
The panel (a) exhibits the locations of the 88 roots for the eigenvalue of maximal magnitude. For the ones of less magnitude the 
typical root distributions are shown on panels (b) and (c) . The panel (d) corresponds to the eigenvalue of minimal magnitude 
where the spirals are converted and the roots located along two ovals. 



The case of -^ = || « ^ is presented in figure 3. In the case of maximal A^ the roots are arranged in two groups 
(figure 3a): straight lines with equiangular separations of 7r/3 and the six points resembling hexagon vertices. Few 
points are deviated from the straight lines, what can be explained by relatively small value of A^. Passing to A^ of 
less magnitude the roots regroup from the straight lines into a circle (figure 3b) which subsequently transforms into 
a hexagon (figure 3c) for the eigenvalue of minimal magnitude. The points left beyond the hexagon can be treated as 
remnants of the straight lines. 

We finally present the case of ^ = || « ^ in figure 4. Interplay between the two groups of roots is again the 
case. One group consists of 8 straight lines with equiangular separation of 7r/4, and the other one resembling certain 
closed contour. Changing A^ we observe that roots rearrange from one group into another and vice versa. For the 
eigenvalue of minimal magnitude (figure 4c) the roots are located basically along the closed contour which appears to 
be an octagon. 

Figure 3 (^ = || w ^) and figure 4 (^ = || w |) demonstrate that the root distributions are of high-accuracy 
symmetric with respect to the rotations by 7r/3 and 7r/4, respectively. These observations provoke to think that in 
the limit i^ — > oo, N ^ oo with -§' ^ ji ^ \' 1' ^^^ approximate symmetry can become exact. In this light the 
limiting case of i^ — >■ oo, N ^ oo with -^ — > finite is of special interest. We presume that in this limit the equation 
(31) may comprise rich symmetry properties. We consider this issue as a matter of separate studies. 



VI. CONCLUSIONS 



We have considered the tight-binding model of spinless particles on a honeycomb lattice in magnetic field. The 
corresponding one-particle Hamiltonian turns out to be expressible in terms of the generators of the quantum group 
Uq{sl2)- We have shown that varying the momentum k the cyclic representation of Uq{sl2) can be continuously 
deformed into the highest weight one. Using the functional representation of Uq{sl2) the Harper equation is refor- 
mulated as a system of two entangled functional equations. Tracing out certain symmetry of the entangled system 
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FIG. 3: Root distributions in the complex plane for ^ = |§ ~ j- There are 89 eigenvalues and the 89 corresponding 
polynomials, each of order of A'^ — 1 = 88. For each eigenvalue we have separate set comprising 88 roots of the corresponding 
polynomial. The panel (a) exhibits the locations of the 88 roots for the eigenvalue of maximal magnitude. The panel (c) 
corresponds to the eigenvalue of minimal magnitude. 
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FIG. 4: Root distributions in the complex plane for ■§- ~ ^ ~ \- There are 89 eigenvalues and the 89 corresponding 
polynomials, each of order of A*' — 1 = 88. For each eigenvalue we have separate set comprising 88 roots of the corresponding 
polynomial. The panel (a) exhibits the locations of the 88 roots for the eigenvalue of maximal magnitude. The panel (c) 
corresponds to the eigenvalue of minimal magnitude. 



we reduce the two equations to a single equation (31) admitting solutions in terms of real polynomials. Bethe-like 
equations determining the roots of polynomials are obtained and the roots are solved out numerically for the different 
eigenvalues A and different magnetic fluxes $. Using graphical plots we have shown that the root distributions are 
highly organized and arranged in regular geometric figures. 

We have stated the quantum group Uq^sh) in terms of q-deformed commutation relations (12) which are compatible 
with g-dcformed co-product structure providing tensor product representations. The representations discussed in 
Sections 3 and 4 correspond to onc-particlc states, while the higher representations constructed via the aforementioned 
g-deformcd product structures will be helpful for describing many-body quantum states. 

The procedure described in Section 4 enables to explicitly construct the complete set of orthogonal one-particle wave 
functions |^fcj) {j = 1, . . . , N) at specific values of momenta set by e+5'==° = iq2 ^ where the functional representation 



of the quantum group becomes reduced to polynomial representation. As pointed out in subsection 3.1 the similar 
reduction takes place at other values of k^ given by e+s'^^'' = ±iq2-3o where jo = 0, 1, . . . , A^ — 1. Consequently 
for each aforementioned value of kx we are able to construct the set of one-particle wave functions which can be 
subsequently employed for calculating variety of matrix elements {^ki.jil ■ • • I^fc2 j2)- I^ particular, calculations of 
the matrix elements of the Coulomb operator are in progress. 
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Appendix A 

Consider first the case oi v — even. We then have 

where 2^ is integer, i.e. \ko\ is the multiple of the A^'th part of the Brillouin zone. We then split the integration area 
in (4) into N horizontal strips of the height -i-^S-. Introducing A^-component vectors (/i = A,B) 

*^(fc) = {c^(fc), c^(fc - feo), . . . , c^(fe - Afco + fco)}, (A2) 

and combining ^a and vj/^ into ^ = {'9 a, *s} wc eventually come to (7). 
In the case oi v = odd we have to write (Al) as 

i^ 47r , . , 

i.e. I fco I is now the multiple of the (2A^)'th part of the Brillouin zone. Correspondingly we split the integration area 
of (4) into 2A^ equal horizontal strips and rewrite (4) as 



H 

Ifbz/{2N) 



|Ht(fc)+H(fc)|dfe, (A4) 



where FBZ/(2A^) denotes any separate strip and 



2N 



y^ c^(fc - nko + ko)cg{k - nfeg + fco)- 



2JV 

+ ^ ^n{k)cl^{k - nko + kQ)cg{k - nko)- 



2N 



^ C„{k)c\{k - nkn)cB{k - nk^ + fco), (A5) 



n=l 



with C„(fc) = 5-"e+^'=°i and Cn(fc) = g-"e+*'="^ 

Wc then decouple the summations into 1 ^ n ^ A^ and A^ + 1 ^ n ^ 2A^ and introduce (fi = A, B) 






X^in{k) ^ -j^c^ik - nko + ^o) H 1= — c^ik - A^fco - nfco + fco). (A6) 



Constructing *p = {x^i, X^2. ■ • ■ . X^jvl and * = {*a, *s} we come to (7). 
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Appendix B 

Discussing the case oi v ~ even in subsection 3.2. we wrote down (20) reflecting the action of the group generators 
on a vector / = (/i, /2, . . . , /at). The same relations are valid also for v = odd only if the aperiodic conditions 
fn+N — (— l)"^^/n are assumed. Provided these conditions are accepted, the functional representation cannot be 
constructed as /„ = /(?") any longer. 

In order to overcome this difficulty we introduce the basis where K and K^^ are diagonal. For ly = odd these basis 
states are given by 



where aj = g^J+i and iV = 2J + 1. 

The actions of group generators on Uj read as 

X+Uj ^ K{e-'^'''q^+^ + e-'^^-^-^-^^j+i, (B2) 

X'uj = K-\e+''"''q-' + e+*'="^g^>,_i, (B3) 

Kuj = -e+^'^'^'^q-'^J-^Uj, (B4) 

K~^Uj = -e-'^-'^q^^+^Uj, (B5) 

where q = -q* = e-<^-^im and n = {-qf^'^- 

Remark that due to A^ = odd we have q^ = 1, i.e. the obstacle related to (7^ = —1 is removed by reformulating the 
action of X in terms of q. Correspondingly, the relations (B.2) - (B.5) require the periodic conditions uj+at = Uj 
which on the other hand are really met by (B.l). Physically this is related to $ — >■ <I>o — $. 

We now pass to another basis Vj set by 

Af 

Wj=A^~5^g^%„. (B6) 

n=l 

Then the relations (B.2) and (B.3) turn into 

K-ijC+Wj = e-*'="ig^i;j„i + e-''"''qh>,+i, (B7) 

kX-Vj = e+''"''q-^-^Vj+i + e+''"'''q^-^Vj_i, (B8) 

which are the analogues of (20a) and (20b). 

Every A^-component column can be expanded as 

M 

: = (/)iui + 021^2 H + i'NVN, (B9) 

In ) 

so the action of group generators on / is now understood as acting on the columns Vj . 

Alternatively, the action of group operators on / can be understood as modifying the corresponding coordinates 
4> = {(pi, . . . , (Jjn)- Provided X^ act on Vj as set by (B.7) and (B.8), the equivalent actions on (p are given by 

K-i(X+0), = e-*'=°igJ"+Vj+i + e-"="^g^-V,-i, (BIO) 

k(X-0)j - e+^'^-^r-'^,-! + e+'^^'^r^V'j+i. (BU) 

Introducing here 

0n = fiKe-i^'y'^r), (B12) 
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we come to the functional representation 

X+f{z) = e~i''-''qzf{qz) + e+^^^^q-^ zf{q~^ z), (B13) 

X-f{z) = e+i'^^'^z-^fiq-'z) + e-^'^^'^z" V(g^). (B14) 

Setting e~5'=='' = i, the eigenvalue equation appears as 

az) + tqzaqz) - iq-^z^(q-^z) = ^C(^), (B15) 

C(z) + iz-\{qz) - tz-\{q-'z) = E^z), (B16) 

and coincides with (30). 

Appendix C 

Performing the replacement z — > — z^^ in (28) we obtain 

f{-z-^) - e+i'^-'^q+h-'^fi-qz-^) - e~i*'-''q~^ z"^ f{q-^ z'^) = X{iz)~^ f{z). (CI) 

Multiplying both sides by \{iz)'^ we pass to 

\{izYf{-z-^) - e+i^'^''q+^ z-^\{izY f{-qz-~^) - e-^^-'^q'^ z-^\{izY f{q-^ z'^) = \^f{z). (C2) 

where the first terms of left hand side can be replaced using (28). The second and third terms can be replaced in 
accord with 

X{izr f{~qz-^) = q-f{q-'z) + g"e+*'^-^-^z/(z) + q^e-i'^^'^q-hfiq-'z) (C3) 

X{iz)'^f{-q-^z-^) = q-'^f{qz) + q-'^e+i''-''q+hf{q^z) + q^'^e-^'^'^'q+hfiz) (C4) 

which can be obtained from (28) hy z ~> qz and z — > q~^z, respectively. 

Performing these manipulations in (C.l) and taking into account g"e"*"'''"=° = —1 we obtain 

(3 - A2)/(z) + €+'''-'' qfiq^z) + e+i''-''{q+h + q-h~^)f{qz) + 

+e-i''-''{q-h + q+h-^)f{q-^z) + e-'''-''q-\f{q-^z) = 0. (C5) 

which coincides with (25a). Reformulating for F{z) = (iz)"/(— z~^) we come to (25b). 
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